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Abstract. We analyze an efficient frequency estimation scheme that is applied 
to measure the unknown frequency of an atomic state in Ramsey spectroscopy. 
The scheme is employing appropriate combinations of uncorrelated probe atoms 
and Greenburgur-Horne-Zeilinger (GHZ) type correlated probe atoms to estimate its 
frequency. The estimation value of frequency is obtained through the Bayesian analysis 
of the final measurement outcomes. The proposed scheme allows us to obtain better 
precision than the scheme without quantum correlation and it also prevents us from 
ambiguity in the frequency estimation procedure with GHZ correlations only. We show 
that the scheme can beat the shot-noise limit and, in addition, it is found that there 
is the trade-off relation between the precision of the frequency estimation and the 
decoherence rate in the atomic states. 
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1. Introduction 

Measurement of a physical system with an arbitrary precision is an important task in 
many fields of experimental physics and technologies. The ultimate goal of the precision 
measurement is to attain the highest accuracy when a physical system is measured at a 
cost of given resources. In a real situation, the aim is needed to be achieved under the 
constraint of limited resources such as the finite number of trials and the limited duration 
time of the experiment. This situation is also true in the case of Ramsey spectroscopy 
experiment which is a scheme to measure the dynamical evolution of an atomic state 
and estimate the resonance frequency of the atomic state [I]. Although precision in the 
measurements can be improved with the number of trials, it is also true that the limited 
resource restricts the possible number of repetitions of the measurement. Therefore, it 
is required to devise a method how to utilize the given resources as efficient as possible. 

A possible method to achieve a higher accuracy in the precision measurement 
is to employ quantum entanglement. The examples are to use coherent-squeezed 
state in gravitational wave detector [2], number entangled states of photons (so- 
called NOON states) in Mach-Zehnder interferometer [3] and GHZ states of atoms in 
frequency estimation with Ramsey spectroscopy (4, 51. In contrast to the conventional 
measurement scheme that allows us to attain the shot-noise limit only, 66 rs-/ i /Vn 0, 
the scheme with entangled states helps us to obtain the Heisenberg limit, 66 ~ 1/N 
Mi where 6 is an unknown parameter to be measured and N is the number of 
employed resources. However, the precision measurement scheme is sometimes restricted 
if a prior knowledge about the value of a parameter to be measured is not provided. 
In the case of phase estimations, there is recent work about an estimation scheme that 
is performed by using entangled coherent states without a prior knowledge about an 
unknown phase |10j . The reason that the precise estimation of a completely unknown 
phase is possible is that without any prior knowledge about the unknown phase, we 
know that phases lie in [0, 27 t) . In general, the precision measurement scheme is only 
useful when we try to achieve a better precision on vaguely known values. If we do not 
have any knowledge on the measurement values, the scheme does not work well. 

Especially, when there is periodicity of probabilities in each measurement, the 
measurement data cannot be used to specify the single value of phase estimation and 
results in more than one estimation value over the period. Thus, if we do not know the 
range of the parameter sharp enough, it is not possible to estimate the precise value of its 
phase and leaves ambiguity. Since the period of probabilities for the measurement with 
entangled states is even shorter than that of unentangled states, entanglement-based 
measurements give rise to a more serious ambiguity problem in the determination of 
an estimation value. It can be shown that in the scenario having a prior knowledge on 
the periodicity, the measurement scheme with highly entangled states merely attains an 
equal precision to that with uncorrelated states in the frequency estimation. Therefore, 
there will be no benefit of using entangled states in the case with limited a priori 
knowledge provided. Here, we propose measurement schemes that allow us to overcome 
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the ambiguity of 2n periodicity and, at the same time, obtain a better precision beyond 
the shot-noise limit using quantum entanglement. 

There have been prior discussions that, when highly entangled states are used, 
periodicity of probabilities in each measurement outcome sometimes prevents us from 
determining a unique estimation value and it can be overcome in idealized cases 
mm |,13j . In the works, two different settings are considered (i) using GHZ states 
with particle numbers 1, 2,4,..., 2 P_1 , and (ii) using combinations of an uncorrelated 
state and GHZ states in order to avoid ambiguity. Each scheme has been discussed 
in interferometric systems for phase estimation and the introduced decoherence model 
was photon losses. In this paper, we consider more general situation with a different 
decoherence effect that occurs in trapped ions. We compare the results of the two 
schemes in the absence of decoherence and in the presence of decoherence, and conclude 
that while the first scheme allows us to achieve the Heisenberg limit in the ideal case, 
the second scheme becomes more advantageous as the decoherence rate increases and 
attains a sub shot-noise precision. 

We organize this article in the following. In section 2, we review frequency 
estimation procedure in Ramsey spectroscopy when the atoms are n n correlated and 
correlated under decoherence effects. Furthermore, we bring up an ambiguity problem 
that occurs in the frequency estimation. In section 3, we analyze two different schemes 
that overcome the ambiguity problem and, at the same time, improve precision by using 
entanglement, and compare the attainable precision of the proposed schemes with that 
of the conventional scheme using no entanglement. Finally, we summarize our works in 
section 4. 

2. Ramsey Spectroscopy 

2.1. Standard frequency estimation scheme under decoherence 

Ramsey spectroscopy is the measurement technique of transition frequency ujq between 
the internal states of two-level atom [T]. The system can only take one of two energy 
levels as they can denote the ground state | g) and the excited state |e). Standard Ramsey 
spectroscopy is operated as following. Initially, prepare N ionized atoms confined in an 
ion trap [IT] and by optical pumping technique [HS], change the state of atoms into the 
same ground state, |-0o) = | g )- Then, the atoms are applied by 7t/2 pulse which leads 
each atom to be in the superposition state of the ground state and the excited state, 
\ifi) = (|<7) + |e))/\/2- After that, (classical) fields with frequency to are applied to the 
atoms for interrogation time t so that the state changes to [ipf) = (| g) + e~ lAt \e))/\/2 in 
a rotating frame, where A = uj 0 — lu denotes the detuning between the atomic transition 
and classical driving field. Finally, the atoms are applied by the second 7t/2 pulse 
which changes the state to |^ 3 ) = sin(Af/2)|g) + cos(Af/2)|e), and the internal state 
of each atom is measured by scattering light and detecting the fluorescence with a 
photomultipier (PMT). We estimate the transition frequency cjq by counting the number 
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Figure 1 . The operation of Ramsey spectroscopy using a single ion is shown. 
The probability of detecting | g) is sin 2 (At/2) and the probability of detecting \e) 
is cos 2 (At/2) in the ideal case. 


of atoms in the excited state |e) and using the probability of detecting the excited state 
|e) in each atom which is given as 


P(e\uj 0 ) = |(e|t/ 3 )| 2 = cos 2 


1 + cos At 
2 


( 1 ) 


The probability of detecting the ground state is given as P(g\u>o) = 1 — P(e|u;o) = 
sin 2 (Af/2). The operation of Ramsey spectroscopy is shown in Fig[T} 

From its functional form, the probability P(e|cuo) can be used to estimate the 
frequency ujq of the atomic state within a certain precision. The statistical fluctuation 
in the estimation which is associated with the given probability P(e|cuo) can be obtained 
by using Cramer-Rao inequality mm, 


59 > 



( 2 ) 


where F(9) = J2i pp|g) the Fisher information and v is the number of 

repetition of trials. P(i\9) is the conditional probability of obtaining a result i with 
a given parameter 9 and the summation is taken over all possible results. The Fisher 
information is a quantity that measures how much information can be obtained when the 
parameter 9 is changed infinitesimally. It is known that the lower bound of the Cramer- 
Rao inequality can be achieved asymptotically by the maximum likelihood estimator 
[18]. The Fisher information plays an important role in estimation theory in that it 
gives a lower bound of statistical fluctuation. 

In Ramsey spectroscopy, the unknown parameter to be estimated is the transition 
frequency u; 0 between two internal states \g) and \e). Using the fact that possible results 
are the excited state |e) with the probability P(e|o;o) and the ground state | g) with the 
probability P(g\u>o), the Cramer-Rao lower bound is easily obtained as 


|<5a;o| = 


(3) 


y/NTt' 

where v = NT ft is the repetition of trials and T is the total experimental time. Thus, 
the precision with N atoms is proportional to l/\ /TV, which is so-called the shot-noise 
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limit [ 6 ]. The shot-noise limit is originated from quantum projection noise. Note 
that the precision |<5u;o| is independent of the transition frequency coq with ignorance 
of decoherence effects. 

So far, we have discussed in the ideal case. However, in a realistic experiment, 
it is necessary to include the effects of decoherence. As a realistic case, we introduce 
a model that was proposed in [19]. In this model, dephasing of individual ions is the 
main decoherence, which is caused from collisions, stray fields, and laser instabilities. 
The effects of decoherence on a single ion p can be described as the master equation in 
Lindblad form [19, 20]: 


p(t) = iA(p\e)(e 


7 

|e)(e|p) + -{a z pa z - p) 


(4) 


where 7 is the decay rate and a z = \g)(g\ — |e)(e| denotes a Pauli spin operator. Eq. ([I]) 
is written in a rotating frame. The effect of decoherence is shown in the broaden signal 
(jl| of a single ion: 


Pj(e |cuo) 


1 + cos(At)e 7t 
2 


(5) 


and again using Cramer-Rao inequality Q, the corresponding precision changes to 


l&M = 


1 — cos 2 (At)e 2lt 
\| NTte~ 2 ^ sin 2 (At) 


( 6 ) 


It can be found that as the decay rate 7 increases, the uncertainty becomes larger. In 
contrast to the ideal case, a precision in the presence of decoherence is dependent to the 
frequency ojq. When At = /c 7 t /2 (k is odd) and t = l/2y < T are satisfied, minimum 
uncertainty is obtained as 


\5oj 


dec 

0 


min 



(7) 


Note that minimum uncertainty is attained when the probabilities of obtaining the 
measurement outcomes | g) and |e) are same. 


2.2. Improvement of precision with GHZ correlation 


We have studied standard Ramsey spectroscopy using N atoms, which allows us to 
obtain a shot-noise limit. It has been proposed that entanglement between atoms 
improves sensitivity of phase. Especially, one of the states that has been proposed 
to improve precision is maximally entangled multipartite state which is so-called 
Greenberger-Horne-Zeilinger (GHZ) state [21] 

| g)® N + |e) 0iV 


I^GHZ) = 


V2 


( 8 ) 


where = \g)\g) ■ ■ • | g) and le)®^ = |e)|e) • • • |e). The GHZ states accumulate N 

times amplified phase information than without entanglement, which results in a better 
precision. Implementation of GHZ states among atoms has been demonstrated by Cirac 
et. al. [ 22 ] and it will be explained in the following. 
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The preparation procedure of GHZ states is, in principle, that after all N atoms 
are prepared in the ground state \g)® N = \g)\g) • ■ ■ | g), the first ion is applied by 
7t/2 pulse to create the state (|g) + |e))|g)|g) • • • |g)/\/2 and the ions are operated by 
a ”controlled-NOT(CNOT)” gate, the first ion as a controlled qubit and the second 
ion as target qubits to entangle the first two atoms, which changes the state to 
( \g)\g) + |e)|e))|g) • • • \g)/\f2. Continuing the operation of CNOT gates, the first ion 
as controlled qubit and other ions as target qubits, the final state becomes the GHZ 
state (\g)\g) ■ ■ • \g) + |e)|e) • • • \e))/\/2 = (| g) 0N + \e)® N )/y/2. Despite theoretical 
straightforwardness of generating GHZ states, it is known that large size GHZ states 
are extremely difficult to create in practice [231 EH 23 E3 2?]- 

After preparation of GHZ states, (classical) fields of frequency to are applied to the 
atoms for interrogation time t , which changes the state to (Ig)®^ + e~ lNAt \e)® N )/\/2. 
Finally, the atoms are disentangled by the second set of controllcd-NOT gates and the 
internal state of the first ion is measured. Again, we estimate the true frequency by the 
measurement outcomes and the probability of detecting |e), 

/'NAt\ 

(9) 


„. | . 1 + cos NAt 9 

P(e\N,u 0 ) = ---= cos 2 [ —— j 


2 V 2 

which denotes the probability of detecting all N atoms in the excited state. The 
probability of detecting all N atoms in the ground state is given as P(g\N,co 0 ) = 
1 — P(e\N,uj 0 ) = sin 2 (lVAf/2). The principle that GHZ states allow us to obtain 
a better precision is originated from the fact that the phase shift of GHZ states in 
interrogation time is amplified N times than without entanglement in the same period, 
which results in the change of the probability of detecting |e). 

As the previous single ion case, using the Cramer-Rao inequality ([2]), we attain a 
VN improved precision than that of uncorrelated N atoms ([3J), 

^ ~~ ^3v¥ ~~ A Ty/Ti' ^ 

where v = T/t is the repetition of trials. Now, uncertainy of frequency is proportional 
to 1/iV, which is referred to as the Heisenberg limit. As the previous case, uncertainty 
does not depend on the frequency in the ideal case. 

Taking into account decoherence, similarly to uncorrelated states, we can easily see 
that the signal (J9]) of GHZ states with N ions becomes 

1 + cos(iV At)e~ Nlt 


P 7 (e|A^,(u 0 ) = 
and uncertainty of frequency becomes 


I&M = 


1 — cos 2 (A^At)e _2Ar7t 
\| N 2 Tte~ 2N * sin 2 (NAt) 


(ii) 


( 12 ) 


The exponential part in the denominator indicates that if N is too large, uncertainty 
|< 5 co 1 q 6C I goes to infinity, which means that a large size of GHZ state is fragile against 
decoherence. Again, uncertainty depends on the frequency and minimum uncertainty 
is achieved when At = kw/2N (k is odd) and t = 1/2 7 A" < T are satisfied. Minimum 
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uncertainty is obtained when the probabilities of obtaining each result | g) and \e) are 
same. As a consequence, minimum uncertainty is same with that of uncorrelated states 
(Jr]) . Therefore, in the presence of decoherence GHZ states do not help attaining a better 
precision than that of uncorrelated states [T9] . 

2.3. Ambiguity of n-period in the frequency estimation 

Let us consider a realistic experiment of frequency estimation with basic concepts 
of Ramsey spectroscopy. Basically, an estimation process is proceeded based on 
measurement outcomes which depend on the parameter to be estimated and the 
probability of obtaining the experimental data. In Ramsey spectroscopy, measurement 
outcomes are consisted of | g) and \e) only, and the probability of detecting the excited 
state |e) is given in ([!]) for a single atom. The probability ([l]) that we use in an estimation 
process has periodicity, which is shown in Fig [2] A difficulty of frequency estimation 
comes from the periodicity because the periodicity prevents us from distinguishing a true 
estimate of frequency among a number of possible estimates. For example, let us suppose 
that the true frequency (Uo me is zero (we set oj = 0 for simplicity i.e. A = o; 0 .), which leads 
to the measurement outcome |e) at all N atoms. Then, our estimate obtained by the 
measurement outcomes is not ooq = 0 but cj 0 = 2irm/t ( m an integer). Our estimation 
result is not unique. In other words, since two different frequencies (o;o)i = 27T m/t and 
( 020)2 = 27t n/t (m, n integers) give the same experimental data, the data does not allow 
us to choose a single estimate of frequency. Hence, even after obtaining measurement 
outcomes, we are still required to determine a single estimation value. Determination 
of a unique estimate can be possible with the help of a prior knowledge about the 
true frequency Uq. If a provided prior knowledge about the true frequency is narrow 
enough to choose a unique estimate, the ambiguity problem is avoided. If a given prior 
knowledge is not enough narrow, we still have ambiguity in the estimation procedure. In 
short, because of periodicity, even when we try infinitely many measurements, it is not 
possible to determine a single estimation value without an appropriate prior knowledge. 
Thus, we are required to know about the true frequency initially with an enough degree 
of accuracy to choose a single estimation value so that we assume that an arbitrary prior 
knowledge about coo is always given in this paper. 

In this situation, let us suppose that we are trying to improve precision. When we 
use N atoms, the ultimate precision is given as (J3|. Under a restriction of the limited 
number of atoms N and a given total duration time T, a possible choice to improve 
preicision is increasing interrogation time t, since uncertainty of the frequency |&Uo| 
decreases as the interrogation time t increases. Hence, we can increase the time t as 
possible in the only restriction t < T where T is sufficiently large. However, as the 
interrogation time t increases, the period 2tt ft of the signal ([TJ) becomes shorter so that 
it causes ambiguity in determining a unique estimation value because a more number 
of frequencies become a possible estimation value as the period becomes shorter. In 
other words, the shorten period makes it harder to distinguish a true frequency from 



Figure 2. Probability distribution of an excited state as a function of wo when t = 1. 
(a) Periodicity of the probability distribution causes ambiguity in determining a single 
peak among the values in the periodic repetitions. If one knows that ui o € [—37 t,37t), 
the estimate of loq can be taken from any values among ±27r and 0. In order to avoid the 
ambiguity, we need to change the interaction time t three times smaller than original. 
The new choice of interaction time allows us to discriminate the peak value of ujq at 
0 from ±27 t. (b) Probability distibution of a single atom (blue) in an excited state 
and that of GHZ states with 3 atoms (purple) are shown. If we know u>o € [— n, w) 
initially, the best choice of interrogation time is t = 1. If we use a single atom, an 
estimation can be uniquely determined as wo = 0. On the other hand, if we use GHZ 
states with 3 atoms, the estimation becomes ambiguous because our estimation value 
is uio = ±27 t/ 3 and 0. To avoid the ambiguity, we need to choose t three times smaller, 
which leads to the same uncertainty <5 wq with the scheme without GHZ correlation. 


other possible estimates. Thus, while the longer interrogation time t allows a better 
preicision, the increased interrogation time t requires a narrower prior knowledge about 
the true frequency coq in order to determine a unique estimation value. 

Let us suppose that we know that the true frequency ojq lies in the interval [0, tt/L) 
where L is a positive number. Under the prior knowledge, the largest possible value of 
interrogation time t that allows us to determine a single peak is L. Then, the possible 
maximum interrogation time t = L leads uncertainty to be 


|<5a;o| = 


(13) 


Vntl 

This is minimum uncertainty with the given prior knowledge loq € [0,7 t/L). It implies 
that minimum uncertainty of frequency estimation depends on the accuracy of a prior 
knowledge about the true frequency. In other words, if we know the true frequency more 
accurately at the beginning, we can obtain a better precision. 

As previously discussed, quantum entanglement between atoms help us to improve 
precision of frequency estimation. Let us consider introducing quantum entanglement 
between atoms by exploiting GHZ states under the same prior knowledge that ojq G 
[0, 7 t/L). Similar to the unentangled state case, the signal of GHZ states with N atoms 
has a period 27 r/Nt from (§■ Therefore, the largest possible t that allows us to choose 
a single estimation value is L/N. Eventually, uncertainty obtained by using GHZ states 
is same with (13). The reason that uncertainty with uncorrelated state and with GHZ 
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state is same is that the period of signal with GHZ states is N times shorter than that of 
nncorrelated state so that we need to use a shorter interrogation time to avoid ambiguity 
in the estimation process. Consequently, in the consideration of a prior knowledge about 
the frequency, exploiting GHZ states does not improve precision, which is shown in Fig 
[2j However, it does not mean that GHZ states are useless in improving the precision 
of frequency because it obviously has a better sensitivity locally. We introduce useful 
schemes that utilize the potential power of GHZ states and improve precision, avoiding 
ambiguity. 


3. Schemes to overcome ambiguity of 7r-period 


In this section, we analyze two measurement schemes to overcome ambiguity in 
frequency estimation of atomic states. Basic settings of the schemes are originally 
proposed by Pezze et. al. mi and Gkortsilas et. al. ra. However, their schemes are 
limited either in the case of ideal measurement or in the case with photon losses. Here, 
we apply the schemes to frequency estimation with atomic GHZ states in the presence 
of decoherence. Decoherence effects on atomic states are different with photon losses 
in photonic states. Since photon losses can be recognized by comparing the number 
of detected photons with that of input photons, one is able to exclude measurement 
outcomes where photon losses occurred. Thus, one takes perfect measurement outcomes 
only in the estimation process. On the other hand, detection of internal states of atoms 
does not give any information whether decoherence had an effect on the atomic states. 
Therefore, all the measurement outcomes should be taken into account including atoms 
that were affected by decoherence. As such, decoherence in atomic states is different 
with photon losses. 

The first scheme is using GHZ states with 1,2, 2 2 ,..., 2 P ” 1 numbers of particles 
in order. It is known that this scheme allows phase estimation with the Heisenberg 
limit in the ideal case. The second scheme is using combinations of an uncorrelated 
state and GHZ states. The two schemes behave differently and give different precisions 
depending on the degree of decoherence effects. We analyze the two schemes in frequency 
estimation, considering decoherence effects on atomic states. 

Before we introduce details of the schemes, let us review Bayesian analysis. Let 
us suppose that we know that the true frequency lies in [0,7 r/L) again, and we choose 
t = L which leads to the best sensitivity avoiding ambiguity for an uncorrelated state. 
According to the Bayes’ theorem, the posterior probability distribution P p (u> 0 \N T , X), 
given total resources Nt and a data set of measurement outcomes X = (aq, # 2 , ■■■, % T ), 
when x n E {g, e}, is given as 


P p (u> 0 \N T ,X) 


P(X\N t ,u 0 )P(u 0 ) 
P(N t ,X) 


(14) 


where P(cno) is a prior distribution of uiq and is completely flat for an unknown 
frequency, P{ojq) = L/ir , and P(X) is treated as a normalization constant. The 
posterior probability P(uj 0 \N t , X) is a conditional probability of frequency u> 0 based 
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Figure 3. Scheme with GHZ states with particle numbers N e = 1,2,4..., 2 P ~ 1 (a) 
Preparation of GHZ states with 1,2 and 4 atoms. The connection line between atoms 
denotes GHZ correlation, (b) Using GHZ states with 1,2 and 4 atoms, the overall 
probability distribution (blue solid line) has only a single peak. WQ rue f = 0. A prior 
knowledge about frequency is ujq £ [—7r,7r). Orange, red and green dotted lines denote 
probability distributions of using GHZ states with 1, 2, and 4 atoms, respectively. 


on a data set A" of measurement outcomes obtained by using Nt resources. We 
choose the maximum of the distribution as estimator uTq and uncertainty |&no| given 
by dcu 0 P p (uj 0 \N T , X) = 0.6827 when P p (lu 0 \N t ,X) goes to normal distribution. 

Bayesian analysis allows us to obtain the estimation value and associated uncertainty 
by using measurement outcomes. 


3.1. GHZ states with particle number 1,2,4, ...2 P 


-1 


Let us consider the first scheme : exploiting GHZ states with particle numbers 
N e = 1, 2,4,..., 2 p_1 . Note that we have assumed that we know that 0 < cUo rue < 7r/L so 
that we have fixed the interrogation time t = L. Now let us suppose (no me = 7t/2L for 
simplicity. In this case, if we use a single atom, the probability of detecting the ground 
state and that of the excited state are same. Thus, the precision is maximum in the 
presence of decoherence. However, if we use GHZ states with even number of atoms, 
the measurement outcomes are either all | g) or all \e) where the precision is minimum 


for the realistic case, which is implied by (12). Thus, for GHZ states with the even 


number of atoms, we need to apply an additional phase shift 7t/2 in order to change the 
overall phase shift into one that gives minimum uncertainty where the probabilities of 
detecting | g) and |e) are same. Then, the effective frequency for GHZ states with the 
even number of atoms is ujq + 7t/2 N e . 

Asymptotically, when we use GHZ states of N atoms (If N = 1, a single atom state) 




11 


with the number of trials v in frequency estimation, we obtain vP(g\N,U)Q S ) number of 
the ground state and uP(e\N, cug®) number of the excited state. Due to the additional 
phase shift, the effective frequency cuq® is (Uq ™ 6 + tt/2N for GHZ states with even N 


atoms and (Uo rue for a single atom. Therefore, by using the Bayes’ theorem (14), when 
we use GHZ states of particle number N e = 1, 2,4,..., 2 P_1 with repetition v = T/L, the 
probability distribution in the ideal case becomes 

P p (u 0 \N t = £2 fc ,X) oc )jP(u 0 \2 k ,ey p W k ^ x P(u 0 \2 k , g y p W\<) 

k =0 k =0 

P - 1 

= [cos 2 ((Uot/2) sin 2 (o;ot/2)] 1 ' //2 [cos 2 {2 k uj^t / 2 — 7t/4) sm 2 {2 k ujot/2 — tt/4)]^ 2 

k =1 


~ e - J/ ( 4 P - 1 )G°*-f ) 2 / 6 


~ e 


-vN 2 (ujot—^) 2 /6 


(15) 


where P(uo\N e ,x ) is dehned by the Bayes’ theorem (14). We have used gaussian 
approximation [sin 2 ((Uoi/2)] m [cos 2 ((Uot/2)] n ~ f or l ar g e m + n, where 

(hot = 2tan _1 (^/m/n) is the maximum point, and Nt = I^k=o ^ = 2 P —1. An advantage 
of using GHZ states with geometrically increasing number of atoms is cancellation of all 
peaks except the central one, which is shown in Fig[3j As a consequence, the probability 
distribution gives uncertainty of frequency which achieves the Heisenberg limit, 


|h(Uo| 


V3 

sj.vNj,L 2 


vg 

NtVtl' 


(16) 


In the realistic case, P(x\N,u> 0 ) are replaced by P y (x\N,oJo). The results are shown in 
Fig [5j Since it achieves the Heisenberg limit in the ideal case, apparently this scheme 
is helpful to utilize GHZ states under a prior knowledge and to improve precision. 
Nevertheless, this scheme is extremely fragile against decoherence because it exploits a 
large size of GHZ states. The result of the first scheme in the presence of decoherence 
indicates that it is required to use a number of uncorrelated atoms that are more robust 
against decoherence than GHZ states. Therefore, we need to introduce another scheme 
that is more robust against decoherence in a realistic situation. 


3.2. Combination of different correlations 

The second scheme is to use appropriate combinations of uncorrelated atoms N u and 
the p copies of GHZ states with N e atoms. Here, the uncorrelated atoms play a role 
in suppressing other periods except one period where the true frequency exists, while 
GHZ states play a role in improving sensivity, which is shown in Fig [4j Similar to the 
previous scheme, the probability distribution with repetition v is asymptotically written 
as 

P p (u 0 \N t = N u + pN e , X) oc 

P((Uo|l, e) t/JVuP ( e l 1 > a; o rue ) x P(u 0 \ hg y N »r(M™) 
xP(u 0 \N e , e yp p ^\ AG< ue ) x P(u 0 \N e , g ypP(9\N B ,^) _ ( 17 ) 
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To minimize the standard deviation of the probability distribution, we need to choose 
optimal N u , N e and p numerically. For numerical optimization, first of all, we obtain the 
probability distribution in ( [I7| by using ((Tj) and ([9j) in ideal case or by using (J5]) and ( [TT| ) 
in realistic case. Then, for a proper value of repetition v and fixed Nt = N u + pN e , by 
changing N u ,N e and p appropriately, we find values of N u ,N e and p that minimize the 
standard deviation of the probability distribution. Then, we iterate the same procedure 
for different values of N t- After the numerical optimization, only one peak survives 
where the true frequency exists. 

Let us consider the same situation with the first scheme, a prior knowledge that 
(UQ me G [0, 7 t/L), fixed t = L and the true frequency u>l Tue = 7 t/2 L. In this case, we 
only use GHZ states with odd number of atoms because GHZ states with odd number 
of atoms attain minimum uncertainty at cjQ rue = 7t/2 L in the presence of decoherence. 
At (UQ me = 7 t/2 L, the probability of detecting | g) and that of detecting \e) are 1/2 
for both uncorrelated states and GHZ states with odd number of atoms. Substituting 
4 rue = 7 t/2 L into (17), in the ideal case, the asymptotic probability distribution with 
the combination of N u uncorrelated atoms and p copies of GHZ states of N e atoms and 
repetition v = T/t can be obtained as 

P p (ui 0 \N T — N u + pN e , X ) 

= [cos 2 (a;oi/2) sin 2 (u Q t / 2)} l/ Nu ^ 2 [cos 2 (N e u 0 t/ 2) sm 2 (N e u 0 t/2 )] up / 2 

~ Ce -( Nu + pN e ' > ( UJ ° t -' n '/ 2 ' )2 / 2 ( 18 ) 


where C is a normalization constant. The gaussian approximation is valid when the 
probability distribution is well-localized around 7 t/ 2 . If N u ,N e and p are numerically 
optimized so that the probability distribution has a single peak, the approximation is 
valid. After the numerical optimization, uncertainty of frequency is obtained as 

|5Wo1 = \I v(N u + pN*)V = N u +pN2)TL (19) 

which is the consistent with the Cramer-Rao bound. It implies that the Cramer-Rao 
bound is saturated by Bayesian analysis if there is no ambiguity in the estimation 
process. Uncertainty obtained by the numerical optimization of N u , N e and p is shown 
in Fig [5} In the realistic case, P(x\N, coo) are replaced by P-y(x\N,co 0 ) in ( |T7] ) . 

In the absence of decoherence, exploiting the GHZ states with 1, 2,4,..., 2 P_1 atoms 
allows us to achieve the Heisenberg limit as known. A drawback of this scheme is 
that as the decoherence rate 7 increases, uncertainty becomes larger as shown in Fig 
[5| Especially, at 7 = 0.1 for large Nt, uncertainty is larger than that of uncorrelated 
states. In other words, since a large size of GHZ state is fragile against decoherence, 
this scheme is not practical when the decoherence rate is large. On the contrary, when 
we use appropriate combinations of an uncorrelated state and GHZ states, since we 
use a number of uncorrelated atoms, we can expect the robustness against decoherence. 
Indeed, Fig [5] shows that the scheme with combinations of an uncorrelated state and 
GHZ states is robust against decoherence. More generally, it is noticed that the best 
scheme depends on decoherence rate 7 . In the ideal case, the first scheme attains 







13 


(a) 




CD 


(c) 



Figure 4. Scheme with combinations of uncorrelated atoms and GHZ states. L = 1, 
t = 1. (a) Preparation of a combination of three uncorrelated atoms and two copies 
of GHZ states with two atoms. The connection line between atoms denotes GHZ 
correlation, (b) An appropriate combination (N u = 7, N e = 3 and p = 4) results in 
a single peak, (c) An inappropriate combination (N u = 1, N e = 3 and p = 2) gives 
rise to a probability on the region that does not include a true value. Green dashed 
line represents a probability distribution from p copies of GHZ states with N e atoms 
and red dotted line represents a probability distribution from uncorrelated N u atoms. 


Heisenberg limit, which is the best sensitivity among onr schemes. As the decoherence 
rate 7 increases, however, the second scheme becomes better. Tims, depending on the 
decoherence rate 7, we need to choose an optimal scheme. 

In the second scheme, we have used combinations of uncorrelated atoms and 
GHZ states with fixed number of atoms for simplicity. However, it is possible to use 
combinations of uncorrelated states and various sizes of GHZ states. Indeed, the first 
scheme is one of the combinations of various GHZ states and the first scheme gives a 
better precision when the decoherence rate 7 is small. Therefore, it can be considered 
to use different various combinations depending on the decoherence rate 7. 

We have assumed that the true frequency u> 0 is 7 t/2 L in the both schemes because 
the best sensitivity is attained at n/2L for uncorrelated atoms in the realistic case. If 
the phase shift u 0 L is not tt/2, we need to modify the phase shift to be tt/2 in order 
to attain minimum uncertainty, which is realizable by a feedback mechanism PEJH3]. 
A feedback mechanism is implemented by the Bayesian analysis on the collected data. 
After each trial, we apply an additional phase shift which is determined by the Bayesian 
analysis to make the total phase shift to be n/2. For every trial, we repeat this process 
and then total phase shift becomes 7 t/ 2 which gives minimum uncertainty of frequency. 

One of the advantages in the second scheme is that it does not require a large size of 
GHZ state which is experimentally difficult to generate and has a short coherence time 
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(c) (d) 




Figure 5. The comparison of precision with various decoherence rates 7 . We fix 
L = 1, t = 1 and Wq ™ 6 = n/2. (a) In the absence of decoherence (ideal case), 

exploiting GHZ states with 1, 2,4, 2 P ~ 1 atoms (orange) achieves the Heisenberg 
limit and the uncorrelated states (red) attain shot-noise limit and combinations of 
uncorrelated atoms and GHZ states (blue) achieves the sub shot-noise limit. (b)(c)(d) 
As the decoherence 7 increases, the scheme with GHZ states with 1 , 2,4,..., 2 P_1 atoms 
losses its advantage since GHZ states with a number of atoms are fragile against 
decoherence. However, the scheme with combinations of uncorrelated states and GHZ 
states beats the shot-noise limit even in the presence of decoherence, (b) 7 = 0 . 01 . (c) 
7 = 0.05. (d) 7 = 0.1. 

[23UMU25U261I27]. Indeed, the largest GHZ states that we have used in the numerical 
optimization are those of 5 atoms in both the ideal case and the realistic case. It 
indicates that our scheme is practical as well as useful. Nevertheless, since generating 
GHZ states requires extremely delicate experiment and the coherence time is very short, 
realization of our scheme may necessitate more advanced experiment devices and skills 
than current. 


4. Conclusion 

Periodicity of probability distribution causes ambiguity in the frequency estimation 
process. Furthermore, in the consideration of a prior knowledge, exploiting GHZ states 
does not help improving precision of frequency estimation when one uses only GHZ 
states with the same number of atoms. In order to avoid the ambiguity and utilize GHZ 
states for improvement of precision, we implement two different schemes. The first 
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scheme is employing GHZ states with 1,2,4,..., 2 P_1 atoms, which improves precision 
significantly in the ideal case. Nevertheless, since a large size of GHZ states is fragile 
against decoherence, the first scheme is no longer advantageous in the realistic case. The 
second scheme that is robust against decoherence is exploiting appropriate combinations 
of uncorrelated atoms and GHZ states. As the decoherence rate 7 increases, it is shown 
that the sensitivity of second scheme can be better than that of the first scheme. In 
addition, we conclude that the scheme with combinations of uncorrelated and GHZ 
states allows us to achieve a sub shot-noise precision in the presence of decoherence. 
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